ANL-HEP-PR-12-xx 



Single- variable asymmetries for measuring the 'Higgs' boson spin 

and CP properties 

Radja Boughezal/'Q Thomas J. LeCompte/'Q and Frank Petriello 1,2, @ 

1 High Energy Physics Division, Argonne National Laboratory, Argonne, IL 60439, USA 
bJQi 2 Department of Physics & Astronomy, 

Northwestern University, Evanston, IL 60208, USA 

Abstract 



J2' 



We introduce a class of asymmetries sensitive to the spin and CP properties of the new boson 
discovered by the ATLAS and CMS experiments. These asymmetries can be measured in the 
four-lepton final state, and are defined by integrating the invariant masses of the lepton pairs over 



specified ranges. We outline a program of measurements using initial LHC data to determine the 
quantum numbers and coupling structure, provide analytic expressions for decay widths in several 
representative models, and discuss what can be determined using the available data. As examples, 
we show how the combination of ATLAS and CMS data already disfavor certain spin-2 couplings, 
and discuss how further data will allow for discrimination of a pure CP-odd scalar from a CP-even 
hypothesis. 
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I. INTRODUCTION 



The recent discovery of a new boson by the ATLAS and CMS collaborations at the 
LHC [l|, 0] has ushered in a new era in particle physics. The future program of the LHC, 
and the next stage of experimental studies in high energy physics, will be largely devoted 
to measuring and understanding the properties of the new state in order to determine the 
underlying theory from which it arises. The initial data provides only a hazy glimpse at the 
properties of the new particle. Observation of its decay into two photons indicates that it 
cannot be a spin-one state, according to the Landau- Yang theorem |3( . Initial measurements 
of its branching ratios into various final states indicates that its couplings are consistent with 
those predicted for the Standard-Model Higgs boson, as determined by the experimental 
collaborations and by several independent analyses |4j-l6j. The slight excess observed over 
Standard Model predictions in the 77 final state has already received explanations both 
within [7] and beyond |8-18| the Standard Model. Significant work will clearly be needed 
to sharpen our picture of the new state. 

Two of the first quantities requiring determinations are the spin and CP properties of 
the observed particle. These can be measured through a variety of ways in multiple final 
states. An initial attempt to determine the CP properties assuming a spin-zero state has 
been made [l9j], and a discussion of how to disentangle the various spin and CP possibilities 
once more data is taken was presented in Ref. 20J. As the largest experimental excesses 
are in the 77 and ZZ — > 41 final states, where either one or both Z-bosons are off-shell, 
initial studies will focus on these two modes. Two broad categories of techniques exist for 
the measurement of the new particle's quantum numbers. Multi-variate methods input all 
of the kinematic information in an event into a likelihood function that can be used to 
exclude hypotheses for the new state's identity. These ideas have been suggested in the 
literature for both analyzing the new particle's properties and assisting in discovery 
and in particular are heavily used in the CMS analysis of the ZZ final state [2'|. In the 
large-time and large-data set limit, such techniques provide the most sensitivity to particle 
properties, since no kinematic information is neglected. Alternatively, single variables that 
provide sensitivity to properties of interest can be studied. Analyses of this type offer the 
advantages of simplicity and clarity over multi-variate approaches, and can be more easily 
implemented to provide answers quickly. They can also indicate which input variables should 
be used to improve the efficiency of multi-variate techniques. 

We introduce here a powerful single-variable measurement in the four-lepton final state 
that discriminates among both the spin and CP possibilities for the new state. It can be 
used when either one or both intermediate Z bosons are off-shell, and is relatively insensitive 
to background contamination. The idea is simple to explain. Let M12 and M34 respectively 
denote the same-flavor lepton pairs with the highest and lowest invariant masses. We sug- 
gestively call the mass of the heavy resonance decaying to the leptons Mh- M34 must satisfy 
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the following inequality: M34 < Mh — Mi 2 . The fall-off of the M34 distribution as this upper 
limit is approached is sensitive to the spin and CP nature of the heavy resonance. Denoting 
the momentum of the 34-pair in the H rest frame as /3, for a pure CP-even spin-zero state 
this distribution decreases linearly in (3. For a pure CP-odd spin-zero state, it falls off as 
/3 3 . For spin- two states, it falls off as either (3, (3 3 , or /3 5 , depending on the couplings of this 
state to spin-one particles. Measurement of this distribution provides a powerful handle on 
the couplings of this new state. 

The sensitivity of the M 34 distribution to the identity of resonances decaying to a pair 
of Z-bosons has been discussed before in the literature j24j . We sharpen and extend this 
observation by defining a class of asymmetries that provide strong discrimination against 
various hypotheses for spin, CP and couplings. These "high-low" asymmetries are defined 
as follows: 

A _ N(M 34 > M cut ) - iV(M 34 < M cut ) 
Mcut N(M U > M cut ) + N(M U < M cut ) ' 1 j 

For a given M 12 , M 34 e [M[ ow ,M H — M 12 ], where M low is determined by the experimental 
cuts. N denotes the number of events in the indicated range of M 34 , and M cut can be chosen 
to provide discrimination between various hypotheses for the f3 n fall-off. The measurement 
of A is advantageous for several reasons. The variables M 12 and M 34 are exactly those 
used in the discovery of the new state, making it easy to perform this analysis quickly with 
the existing data. The background in this channel is fairly low. With very few events the 
allowed spin and CP combinations are already reduced by measuring A. In this paper we 
investigate this idea in the following ways 

• We provide the decay widths differential in M 12 and M 34 in several reference models 
for use in experimental studies. 

• We outline a program of analyses that can exclude various spin, CP and coupling 
hypotheses. 

• We demonstrate using the public ATLAS and CMS data that several spin-2 coupling 
possibilities are already disfavored. 

A word of caution on the last item is necessary. Even in an an analysis this simple, there 
are details, subtleties and systematic effects in interpreting the data that can only be prop- 
erly treated by the experimental collaborations. Conclusive evidence for or exclusion of a 
hypothesis can only be provided by ATLAS and CMS. Our goal here is only to demonstrate 
the power of the A analysis using the available data as an example. We encourage the 
experimental collaborations to perform a thorough and conclusive study. 

Our paper is organized as follows. We present the differential decay widths for several 
representative spin and CP combinations in Section [Til We discuss the construction of the 
asymmetry in Section lITlj and show how the current ATLAS data can be used to discriminate 
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among the various hypotheses. We discuss our results and outline a program of experimental 
analysis in Section HVl 



II. DERIVATION OF THE DECAY WIDTHS 

We provide here for completeness a derivation of the decay widths for the heavy-resonance 
decay into four leptons for several representative examples. We study the mode H(pjj) — > 
Z{j> ee )Z{j>„) — > e + (p e+ )e~(p e _)/i + (p^ + ) / u~(j^_), and extend the result to the case of four 
same-flavor leptons at the end. The intermediate Z bosons may be either on-shell or off- 
shell. We begin by considering the spin-zero hypothesis for H. The HZZ Feynman rule can 
be written in the generic form iT^ u {p ee) p^). It is straightforward to derive the following 
expression for the differential decay width: 
d 2 T G F M\ 



ds ee ds„ 36n 4 M H 



BW(s ee )BW(s w ) / dPS 2 (pee,P^)T llu (p ee ,p flfl )Tj )r7 (p ee ,p flfl ) 



SeeS,, { g» P ~ ^4 \lT ~ \ ■ (2) 



SflfM 



We have set s ee = p ee = (p e+ + p e J) and = p = ip,+ +P,~) ■ These quantities become 
the M12 and M34 introduced in Section [J when a hierarchy between them is specified. We 
have introduced the following notation for the Breit-Wigner distribution of the Z-bosons 
and the two-particle phase space: 

1 



BW(s) 



(*-Ml)2 + r|Mf 

dPS 2 (p ee! P w ) = J d 4 p ee d 4 p^5(p 2 ee - s ee )5(p 2 ^- s^)5 {4 \p H -p ee -p^), (3) 

where T z denotes the Z-boson width. 

To proceed further, we must specify the form of T MV . We will consider the two cases of a 
pure CP-even scalar and a pure CP-odd scalar. It is possible to also consider a mixed state, 
but we do not pursue that option here. The forms for the interaction vertices in each case 
are given below: 

++ : T^(p ee ,p^) = a 1 M z g^, 

0-+ : T^(p ee ,p^) = a 2 e^Pee, P P^,a/A. (4) 

A has dimensions of energy, while ai and a 2 are dimensionless constants. We note that for a 

Standard-Model Higgs boson at tree-level, a± = g/cw and a 2 = 0. Since we will later utilize 

a shape analysis, the values of these parameters do not matter. Using these expressions it 

is straightforward to derive the following differential decay widths: 

d 2 T(0 ++ ) , in GpM% , v , , a f (3 2 M% 
' W W^CTT BW(s ee )BW(s llfl )(3{3s ee s^+ l 11 



ds ee ds m ™ 72n*M H ™ ™" " 4 
d 2 r((T+) , , 2 G pM^Mft _ TTr/ ,^ T/ , „, 
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We have introduced the relative momentum of the off-shell Z* bosons in the H rest frame, 
in units of M#: 



Mfj + s 2 ee + s 2 - 2M 2 H s ee - 2M 2 H s^ - 2s ee s Ml 
W H ' < 6 > 

We note that (3 = when ^J~s^ = Mh — y/s^- As claimed earlier, the ++ distribution 
decreases linearly as this upper limit is approached, while the h distribution falls off as 

We next derive the decay width for a spin-two state. We write its coupling to Z-bosons 
as iT^'^lpee^nn). It is straightforward to derive the following differential decay width: 

d 2 T G M 2 f 

dg ds = 36 ^m h BW ( g ee) BW(s ff ) J dPS(p ee ,p MM )T a ^(p ee ,p w )T 7 t 5p(r (p ee ,p At/ ,) 

x s ee s^Jg» p -?^\{g u ° -^^\B a ^ 5 . (7) 

B denotes the numerator for a massive spin-two particle propagator. It can be obtained from 



Refs. 25|, |26|]. We must again specify a form of the spin- two particle coupling to Z-bosons. 
A natural choice is to consider a spin- two Kaluza-Klein graviton, which we study first. 
Another approach is to introduce a model-independent parameterization of the coupling, as 
was followed in Ref. j2l|. We will also examine that possibility. 

We denote the spin-two graviton with the notation 2 G . The Feynman rules for this 
state can be obtained from Ref. 
expression: 

d 2 T(2 G ) G f M 2 k 2 



261 ] . The differential decay width is given by the following 



ds ee ds w 288vr 3 M H 



BW(s ee ) BW(s (8) 



k is a constant with dimensions of inverse mass. Its numerical value will not be important in 
our analysis. The full expression for the function Q is lengthy and not especially illuminating. 
We present it for completeness in an Appendix. However, we note here that it does not vanish 
in the (3 — > limit. We demonstrate this by setting J~s^ = Mh — a/s^I, which is equivalent 
to setting /3 = 0. We find 



d 2 T(2 G ) ^ G f M 2 k 2 



ds ee ds m 288tt 3 M h 



BWWBW(g««v^ {40 {M H ^T e -s ee + Ml) 2 }. (9) 



The spin- two graviton decay width therefore decreases linearly with /3 as — > 0. 

A spin-two resonance does not necessarily have to be a Kaluza-Klein graviton. It could 
instead arise as a composite from a strongly-coupled theory, and have different couplings to 
gauge bosons than the graviton. A useful parameterization of spin-two couplings to vector 
bosons was presented in Ref. j2l|. The generic vertex has ten couplings. Of these ten 
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couplings, only three lead to a linear decrease in the M34 distribution as ft — > 0. The others 
lead to distributions that fall off as either ft 3 or ft 5 as ft — > 0. This behavior is simple to 
understand; the coupling that produces a ft 3 behavior is that of a CP-odd "pseudo-tensor" , 
while those which produce ft 5 behavior have additional momentum insertions and therefore 
are produced in a higher partial-wave state. As we have already introduced an example 
with ft 3 dependence above, we choose an effective operator which leads to a ft 5 behavior for 
demonstration purposes: 

r^( Pee , Pw ) = ^(pee-PwY(Pee-PwY {Pee • P^ ~ feeV^ } ■ (10) 



This corresponds to the coupling g\ in Ref . [21[ . We denote this spin-two possibility as 2 



It is straightforward to derive the differential decay width of this state: 

This features a much sharper decrease as ft — > 0. 

We now comment on the extension to the eeee and /x/x/x/x final states. These contain both 
an additional symmetry factor of 1/4 from the identical particles in the final state, and also 
arise from two Feynman diagrams rather than a single one as for the ee/x/x mode. As we are 
interested in a shape analysis, the symmetry factor does not affect our results. It can be 
simply included in the formulae above if needed. Similarly, each of the individual diagrams 
squared gives expressions for the differential decay width identical to those presented above, 
and simply contribute a multiplicity factor of 2, to which our analysis is again not sensitive. 
We now assume that one of the intermediate Z-bosons is close to its mass-shell. This can be 
checked to be true for the majority of the ATLAS data Then, the squared diagrams will 
have one of the Breit-Wigner functions becoming large: BW ~ 1/T 2 . The interference term 
is not enhanced in this kinematic region, and can be neglected. Therefore, the differential 
decay width presented above can be safely used to perform shape analyses in all leptonic 
final states. It is simple to extend the formulae above to account for the interference if 
necessary. 

Before proceeding, we note that this analysis can be extended to higher spin states than 



spin-two. As discussed in Ref. [24j, these fall off quickly as (3 2J ~ 3 , where J denotes the spin of 
the state. We also note that we have not discussed the spin correlations between production 
and decay that occur for the spin-two state. Since we are only measuring invariant mass 
distributions of final-state particles, these will not affect our study. 



III. ANALYSIS 

We begin our discussion of how to test spin and CP hypotheses for the new boson by 
presenting the normalized differential decay widths in Fig. [TJ We have set Mh = 126 GeV, 
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20 25 30 35 

M 34 (GeV) 

FIG. I: The normalized decay width l/T tot x d 2 T / (dM^dMfJ function of M34 for the spin-0 
CP-even and CP-odd states ++ and '", the spin-2 graviton 2 G , and the alternate spin- two state 
2 A . The parameter choices are described in the text. 

and for simplicity have set M 12 = Mz- The lower invariant mass M 34 is shown in the range 
M 34 G [17.5 GeV, M H — M Z ], where the upper limit corresponds to the kinematically allowed 
maximum, and the lower limit to that set by ATLAS. CMS sets their lower limit to 12 GeV. 
The normalization of the decay widths is performed by dividing the expressions presented in 
the previous section by their integrals over the allowed range in M 34 . The shape difference 
induced by the differing j3 n behavior of the various states is clear in the plot. The _+ and 
2 A states are peaked toward lower M 34 , the ++ state is flat, and the 2 G is peaked toward 
high M 34 . 

It is now simple to construct a single variable capable of discriminating between various 
hypotheses. We begin by considering the discrimination between the ++ and 2 A states. 
The two distributions cross at an invariant mass M34 w 26 GeV. It is intuitively clear that 
we should compare the integrals of the two distributions above and below the crossing point. 
We form the asymmetry denned in Eq. ([1]) with M cut = 26 GeV. Choosing M cut to be the 
crossover point can be shown to maximize the sensitivity to the different hypotheses. We 
compute this for the two states under consideration and find the following results: 

AZ 9 (0 ++ ) = -0.078, A s j 6 9 (2 A ) = -0.57. (12) 
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This result for the signal asymmetry must be combined with the background asymmetry 
in order to obtain a prediction that can be compared to the result measured by the ex- 
perimental collaborations. To compute the background asymmetry, we first note that the 
composition of signal and background in the kinematic region of interest is expected to be 
roughly 1:1 from experimental studies [l|, 0]. We further note that it is dominated by contin- 
uum production of ZZ* . We estimate the background asymmetry by running a leading-order 



Madgraph |27| simulation of the pp —> Zl + l~ process with basic acceptance cuts reproducing 
those in the ATLAS analysis. We find an asymmetry of A —0.045. Assuming that 

50% of the measured rate comes from background events, we find the following predictions 
for the combined signal plus background asymmetries: 

A s jg +back (0 ++ ) = -0.060, A S 2 6 9+back (2 A ) = -0.31. (13) 

We can now use the available LHC data to estimate the experimental value of ^4.26 together 
with the statistical error on the measurement. We keep events that are ±5 GeV from the Z- 
peak in order to use the simple choice of M cut = 26 GeV; we discuss this point further below. 
We caution that this is an estimate only. Corrections for detector acceptance, experimental 
resolutions on the invariant masses, systematic errors and other effects must be considered. 
Similarly, events with M 12 ^ Mz should be included to increase the sensitivity to different 
hypotheses. For simplicity we neglect these events, and comment on extending the analysis 
to include them later in the text. Only the experimental collaborations can authoritatively 
perform this analysis. With these caveats mentioned, ATLAS reports 4 events in the higher 
mass bin and 4 in the lower mass bin; CMS reports 1 and 1. for a combined raw asymmetry 
of 0.0 with a statistical uncertainty of approximately 0.28. The 2 A hypothesis is disfavored 
by the data at about the one standard-deviation level. We note that other (3 5 couplings lead 
to similar results; for example, using the effective coupling corresponding to gf^ in Ref. [21I ] 
leads to an asymmetry A%s(g? ) = —0.69, in even more disagreement with the LHC results. 
The message we wish to convey is that decay widths behaving as /3 n with n > 5 are unlikely 
to be compatible with the current data. 

We proceed to consider whether the state h is consistent with the current data. We 
form the asymmetry *4 28 , which provides the strongest discrimination between the ++ and 
0~ + states. We find the following results using the expressions in Section [Til 

A% 8 9 (0 ++ ) = -0.33, ^28 9 (0~ + ) = -0.58. (14) 

Treating the background as before, these become 

A s j 8 9+back {0 ++ ) = -0.31, A s j 8 9+bach {0- + ) = -0.44. (15) 

The combined experimental asymmetry from ATLAS and CMS is —0.40 ± 0.27. At this 
point, both possibilities are consistent with the data. 
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Finally, we compare the 2 distribution with the ++ hypothesis. The crossover point 
for these two distributions is roughly 25 GeV. We compute the A25 asymmetries for both 
states and find the following results: 

A s j 5 g (0 ++ ) = 0.05, A s j£(2 G ) = 0.18. (16) 

Treating the background as before, these become 

A s j 5 9+back (0 ++ ) = 0.06, A s H +back {2 G ) = 0.13. (17) 

The uncertainty on the experimental result is ±0.29, about four times the separation be- 
tween the two hypotheses, so there is currently no discrimination possible between the two 
hypotheses. 

In a more complete analysis, events with M 12 7^ M z would be used. This would double 
the available amount of data. For M12 < Mz, the available phase space for M34 increases, 
indicating that the crossover point also changes. A full analysis should consider M cut as a 
function of M\%. Also, we have selected M cut based on the properties of the signal cross 
section only. It is possible that the chosen values could shift slightly if the combined signal 
and background prediction were used to select M cut . We leave these extensions to a more 
complete experimental analysis utilizing a larger dataset. 



IV. DISCUSSION AND CONCLUSIONS 



We have demonstrated that a simple asymmetry which can be easily formed from the 
available LHC data provides discrimination between the possible spin and CP hypotheses for 
the new boson recently discovered by ATLAS and CMS. The idea behind the measurement 
is simple: as the lesser lepton-pair invariant mass reaches its upper kinematic limit in the 
four-lepton final state, the decay width falls off with a characteristic power of /3. The "high- 
low" asymmetry defined in Eq. (JT]) captures this behavior in a single variable. We have 
shown that even the few events measured in the four-lepton final state are beginning to 
restrict the possible couplings of a spin-two state. An additional benefit to this approach is 
that it facilitates a simple combination of data from multiple experiments. 

These observations suggest the following program of analysis for the experimental collab- 
orations. 

• The experimental collaborations should confirm our expectation that the possible spin- 
two couplings cannot be composed solely of those containing a /3 5 dependence. We 
have demonstrated this by studying two examples from the general parameterization 
presented in Ref . [2l| . A more thorough analysis of the possible coupling structures 
should be performed. 
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• We have focused on events with M 12 ~ M z for simplicity of presentation. A more 
complete analysis would use all of the available data, and would allow M cut to vary 
for events with different Myi- This would increase the sensitivity of the method, and 
should be investigated. 

• The currently available data cannot yet test the CP-odd spin-zero hypothesis. This 
analysis should be pursued as the experimental collaborations observe additional 
events in the four-lepton channel. 

• It is not possible with the current data to distinguish between the spin-two graviton 
and the ++ hypotheses. As more data is collected, this analysis should be revisited. 
Angular distributions will also help discriminate between these hypotheses. 

• With more data it will become possible to extend this analysis to determine whether 
the observed state has admixtures of the various coupling structures; for example, the 
collaborations could consider whether a spin-zero state possesses both CP-even and 
CP-odd couplings. 

In summary, we are excited about the possibility of our proposed simple measurement of the 
properties of the new particle discovered at the ATLAS and CMS experiments. We encourage 
the experimental collaborations to analyze their data using the techniques proposed here. 
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APPENDIX 

We present here the function Q that appears in the differential decay width for a spin-two 
graviton to Z-bosons: 
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j 1 (M z 4 + 12 s ee s^)M H 4 

1 3 s ee s^^ 

2 (3 M z 4 s w - 9 s ee s^ 2 + 3 M z ^s ee + 40 s ee s^ M z 2 - 9 s ee 2 Sflfl ) M H 2 

3 S e e 

2 -42 s ee s w M z 4 - See V s + 7 M z A s^ 2 - s ee 3 s^ + 7 M/s ee 2 

3 s ee s^^i 

2 -34 s ee 2 gwt 2 + 20 See 2 s w M z 2 + 20 s ee s^ 2 M z 2 

3 See 

2 (g ee - s^f (s ee 2 s fMIM + 3 M z 4 5 ee + s ee s^ 2 - 20 s ee s w M^ 2 + 3 Mg 4 ^) 

3 s ee s w M F 2 

1 (S ee ~ S wt ) 4 (2 See S W + M^ 4 ) 
3 See Snn M H 4 




[1] G. Aad et al. [ATLAS Collaboration], [arXiv: 1207.72141 [hep-ex]]. 

[2] S. Chatrchyan et al. [CMS Collaboration], |arXiv:1207.7235l [hep-ex]]. 

[3] L. D. Landau, Dokl. Akad. Nawk., USSR 60, 207 (1948); C. N. Yang, Phys. Rev. 77, 242 
(1950). 

[4] D. Carmi, A. Falkowski, E. Kuflik, T. Volansky and J. Zuoan. larXiv:1207.1718l [hep-ph], 
[5] T. Plehn and M. Rauch, larXiv:1207.6l"08l [hep-ph]. 
[6] A. Diouadi. larXiv: 1208.34361 [hep-ph]. 

[7] J. Baglio, A. Djouadi and R. M. Codbole. la7Xiv:1207.145T1 [hep-ph], 

[8] M. Carena, I. Low and C. E. M. Wagner, JHEP 1208, 060 (2012) [arXiv:1206.1082l [hep-ph]]. 

[9] A. Joglekar, P. Schwaller and C. E. M. Wagner. larXiv: 1207.42351 [hep-ph], 
[10] G. F. Giudice, P. Paradisi and A. Strumia. larXiv:1207.6393l [hep-ph]. 
[11] A. Delgado, G. Nardini and M. Quiros. larXiv: 1207.65961 [hep-ph]. 
[12] M. Hashimoto and V. A. Miransky, larXiv: 1208.13051 [hep-ph] . 
[13] K. Howe and P. Saraswat. larXiv: 1208. 15421 [hep-ph]. 

[14] T. Li, J. A. Maxin, D. V. Nanopoulos and J. W. Walker, l arXiv: 1208. 19991 [hep-ph]. 
[15] R. Sato, K. Tobioka and N. Yokozaki, larXiv: 1208.26301 [hep-ph]. 
[16] Z. Kang, T. Li, J. Li and Y. Liu. larXiv:1208.2673l [hep-ph], 

[17] K. J. Bae, K. Choi, E. J. Chun, S. H. Im, C. B. Park and C. S. Shin. larXiv:1208.2555l [hep-ph], 
[18] H. Davoudiasl, H. -S. Lee and W. J. Marciano. larXiv: 1208 .29731 [hep-ph]. 
[19] B. Coleppa, K. Kumar and H. E. Logan. larXiv:1208.2692l [hep-ph], 

[20] S. Bolognesi, Y. Gao, A. V. Gritsan, K. Melnikov, M. Schulze, N. V. Tran, and A. Whitbeck, 

larXiv: 1208.40181 [hep-ph]. 
[21] Y. Gao, A. V. Gritsan, Z. Guo, K. Melnikov, M. Schulze and N. V. Tran, Phys. Rev. D 81, 

075022 (2010) garXiv:1001. 3396 [hep-ph]]. 
[22] A. De Rujula, J. Lykken, M. Pierini, C. Rogan and M. Spiropulu, Phys. Rev. D 82, 013003 

(2010) iarXiv:1001. 5300 [hep-ph]]. 



11 



[23] J. S. Gainer, K. Kumar, I. Low and R. Vega-Morales, JHEP 1111, 027 (2011) |arXiv:1108. 2274 
[hep-ph]]. 

[24] S. Y. Choi, D. J. Miller, M. M. Muhlleitner and P. M. Zerwas, Phys. Lett. B 553, 61 (2003) 
|hep-ph/0210077| . 

[25] G. F. Giudice, R. Rattazzi and J. D. Wells, Nucl. Phys. B 544, 3 (1999) [hep-ph/9811291 



[26] T. Han, J. D. Lykken and R. -J. Zhang, Phys. Rev. D 59, 105006 (1999) [hep-ph/9811350] . 
[27] J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer and T. Stelzer, JHEP 1106, 128 (2011) 

|arXiv:1106.0522l [hep-ph]]. 
[28] https : //twiki . cern . ch/twiki/bin/view/CMSPublic/Higl2016TWiki 



12 



